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Abstract 

We present a new method for computing complete one-loop amplitudes, including their rational 
parts, in non-supersymmetric gauge theory. This method merges the unitarity method with on- 
shell recursion relations. It systematizes a unitarity-factorization bootstrap approach previously 
applied by the authors to the one-loop amplitudes required for next-to-leading order QCD correc- 
tions to the processes e+e^ ^ Z, 7* ^ 4 jets and pp ^ W + 2 jets. Wc illustrate the method 
by reproducing the one-loop color-ordered five-gluon helicity amplitudes in QCD that interfere 
with the tree amplitude, namely A5;i(l~, 2~, 3+, 4+, 5+) and A5;i(l~, 2+, 3~, 4+, 5+). Then we de- 
scribe the construction of the six- and seven-gluon amplitudes with two adjacent negative- helicity 
gluons, ^6;i(l~,2~",3+,4+,5+,6+) and A7;i(l", 2", 3+, 4+, 5+, 6+, 7+), which uses the previously- 
computed logarithmic parts of the amplitudes as input. We present a compact expression for the 
six-gluon amplitude. No loop integrals are required to obtain the rational parts. 

PACS numbers: ll.15.Bt, 11.25.Db, 11.25.Tq, 11.55.Bq, 12.38.Bx 
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I. INTRODUCTION 



The approaching dawn of the experimental program at CERN's Large Hadron CoUider 
calls for theoretical support in a number of areas. A key ingredient in the quest to find 
and understand the new physics at the TeV scale will be our ability to deliver precise 
predictions for a variety of observable processes. Fulfilling this demand will depend in turn 
on having versatile tools for calculating multi-particle, loop-level scattering amplitudes in 
the component gauge theories of the Standard Model. Tree-level amplitudes provide a first 
but insufficient step. The size and scale- variation of the strong coupling constant imply that 
even for a basic quantitative understanding, one must also include the one-loop amplitudes 
which enter into next-to- leading order corrections to cross sections [1]. An important class 
of computations are of perturbative QCD and QCD-associated processes. Extending the set 
of available processes to -I- multi-jet production, and beyond, will demand computations 
of new one-loop amplitudes in perturbative QCD. 

In this paper we will describe a new approach to computing complete one-loop scattering 
amplitudes in non-sup ersymmetric theories such as QCD. This approach systematizes a 
unitarity-factorization bootstrap approach applied by the authors to the computation of the 
one-loop scattering amplitudes needed for Z ^ 4 jets and pp ^ W + 2 jets at next-to- 
leading order in the QCD coupling [2]. As in that paper, the cut-containing logarithmic 
and polylogarithmic terms are computed using the unitarity method [3, 4, 5, 6, 7, 8] and 
four-dimensional tree-level amplitudes as input. The remaining rational- function pieces are 
computed via a factorization bootstrap, in the form of an on-shell recurrence relation [9, 10, 
11, 12]. (In ref. [2] the rational functions were constructed as ansatze with the assistance of 
the factorization limits, and verified by numerical comparison to a direct Feynman diagram 
computation.) 

The unitarity method has proven to be an effective means of computing the logarithmic 
and polylogarithmic terms in gauge theory amplitudes at one and two loops. In massless 
sup ersymmetric theories the complete one-loop amplitudes may be determined from the four- 
dimensional cuts [4]. This method has been applied in a variety of amplitude calculations 
in QCD [2, 13, 14, 15, 16, 17] and in supersymmetric gauge theories [3, 4, 18, 19, 20]. A 
recent improvement to the unitarity method [7] uses complex momenta within generalized 
unitarity [2, 16, 19], and allows a simple determination of box integral coefficients. (The 
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name 'generalized unitarity', as applied to amplitudes for massive particles, can be traced 
back to ref. [21].) The unitarity method has spawned a number of related techniques, 
include the very beautiful application of maximally-helicity-violating (MHV) vertices to 
loop calculations [8, 22] and the use [23, 24] of the holomorphic anomaly [25] to evaluate the 
cuts. The unitarity method can also be used to determine complete amplitudes, including 
all rational pieces [5, 13, 14, 17] by applying full D-dimensional unitarity, where D = 4 — 2e 
is the parameter of dimensional regularization [26] . This approach requires the computation 
of tree amplitudes where at least two of the momenta are in D dimensions. For one-loop 
amplitudes containing only external gluons, these tree amplitudes can be interpreted as four- 
dimensional amplitudes but with massive scalars. Recent work has used on-shell recursive 
techniques [9, 10] to extend the number of known massive-scalar amplitudes [27]. At present, 
the £)-dimensional unitarity approach has been applied to all n-gluon amplitudes with n — 
4 [17] and to special helicity configurations with n up to 6 [13, 17]. 

The somewhat greater complexity of the /^-dimensional cuts suggests that it is worthwhile 
to explore other methods of obtaining the rational terms. We have additional information 
about these terms, after all, beyond the knowledge that their L>-dimensional cuts are D- 
dimensional tree amplitudes. Because we know a priori the factorization properties of the 
complete one- loop amplitude [3, 28], we also know the factorization properties of the pure 
rational terms. It would be good to bring this information to bear on the problem. This idea 
was behind the 'bootstrap' approach used in ref. [2]. The idea was used to produce compact 
expressions for the Z qqgg amplitudes. However, it was not presented in a systematic 
form, and indeed, for sufficiently complicated amplitudes it can be difficult to find ansatze 
with the proper factorization properties. This shortcoming has prevented wider application 
of these ideas. 

Recent progress in calculations of gauge-theory amplitudes has led us to re-examine the 
bootstrap approach. This progress has been stimulated by Witten's proposal of a weak- 
weak duality between J\f = 4 supersymmetric gauge theory and the topological open-string 
B model in twistor space [29]. (The roots of the duality lie in Nair's description [30] of 
the simplest gauge theory amplitudes.) Witten also made the beautiful conjecture that the 
amplitudes are supported on a set of algebraic curves in twistor space. The underlying 
twistor structure of gauge theories, as revealed by further investigation [23, 31, 32, 33, 34, 
35, 36], has turned out to be even simpler than originally conjectured. (For a recent review. 
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see ref. [37].) The underlying twistor structure was made manifest by Cachazo, Svrcek and 
Witten [38], in a new set of diagrammatic rules for computing all tree- level amplitudes, 
which use MHV amplitudes as vertices. These MHV rules led to further progress in the 
computation of tree- level [9, 10, 27, 38, 39, 40, 41, 42] amplitudes. Brandhuber, Spence, 
and Travaglini [22] provided the hnk between loop computations using MHV vertices and 
those done in the unitarity-based method. This development in turn opened the way for 
further computations and insight at one loop [7, 8, 19, 20, 24, 43]. The remarkable conclusion 
of all these studies is that gauge theory amplitudes, especially in supersymmetric theories, 
are much simpler than had been anticipated, even in light of known, simple, results. Several 
groups have also studied multi-loop amplitudes, and have found evidence for remarkable 
simplicity, at least for maximal supersymmetry [18]. 

Recently, Britto, Cachazo and Feng wrote down [9] a new set of tree-level recursion re- 
lations. Recursion relations have long been used in QCD [44, 45], and are an elegant and 
efficient means for computing tree-level amplitudes. The new recursion relations differ in 
that they employ only on-shell amplitudes (at complex values of the external momenta). 
These relations were stimulated by the compact forms of seven- and higher-point tree ampli- 
tudes [19, 20, 41] that emerged from studying infrared consistency equations [46] for one-loop 
amplitudes. A simple and elegant proof of the relation using special complex continuations 
of the external momenta has been given by Britto, Cachazo, Feng and Witten [10]. Its appli- 
cation yields compact expressions for tree amplitudes in gravity as well as gauge theory [42] , 
and extends to massive theories as well [27]. 

In principle, recursion relations of this type could provide a systematic way to carry out 
the factorization bootstrap at one loop. One must however confront a number of subtleties 
in attempting to extend them from tree to loop level. The most obvious problem is that the 
proof of the tree-level recursion relations relies on the amplitudes having only simple poles; 
loop amplitudes in general have branch cuts. Moreover, the factorization properties of loop 
amplitudes evaluated at complex momenta are not fully understood; unlike the case of real 
momenta, there are no theorems specifying these properties. Indeed, there are double pole 
and 'unreal' pole contributions that must be taken into account [11, 12]. 

In a pair of previous papers [11, 12] we have applied on-shell recursion relations to the 
study of finite one- loop amplitudes in QCD. These helicity amplitudes vanish at tree level. 
Accordingly, the one-loop amplitudes are finite, and possessing no four-dimensional cuts, are 
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purely rational functions. Through careful choices of shift variables and studies of known 
amplitudes, wc found appropriate double and unreal pole contributions for the recursion 
relations, and used them to recompute known gluon amplitudes, and to compute fermionic 
ones for the first time. 

While we will not give a derivation of complex factorization in the present paper, it is 
heartening that no new subtleties of this sort arise in the amplitudes studied here, beyond 
those studied in refs. [11, 12]. The systematization we shall present suggests that a proper 
and general derivation of the complex factorization behavior should indeed be possible. 

In this paper, we focus on the issue of setting up on-shell recursion relations in the presence 
of branch cuts. We describe a new method for merging the unitarity technique with the on- 
shell recursion procedure. As mentioned above, we follow the procedure introduced in ref . [2] , 
determining the cut-containing logarithms and polylogarithms via the unitarity method, and 
then determining the rational functions via a factorization bootstrap. We derive on-shell 
recursion relations for accomplishing the bootstrap. In general, both the rational functions 
and cut pieces have spurious singularities which cancel against each other. These spurious 
singularities would interfere with the recursion because their factorization properties are not 
universal. We solve this problem by using functions which are manifestly free of the spurious 
singularities, at the price of adding some rational functions to the cut parts. These added 
rational functions have an overlap with the on-shell recursion. To handle this situation, we 
derive a recursion relation which accounts for these overlap terms. 

To illustrate our bootstrap method we recompute the rational-function parts of the 
known [47] five-gluon amplitudes. We present all the intermediate steps determining the 
rational functions of one of the five-gluon amplitudes, in order to underline the algebraic 
simplicity of the procedure. As a demonstration of its utility, we also compute two new 
results, the six- and seven-gluon amplitudes with two color-adjacent negative helicities. We 
present the complete six-gluon amplitude in a compact form. These results have all the 
required factorization properties in real momenta, a highly non-trivial consistency check. A 
computation based purely on the unitarity method, that is to say based on full I?-dimensional 
unitarity, would provide a further check. 

This paper is organized as follows. In the next section, we review our notation and the 
elements entering into a decomposition of QCD amplitudes at tree level and one loop. In 
section III, we derive a new on-shell recursion-based formula for general one-loop ampli- 
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tudes. In section IV, we review the relevant known amplitudes, and pieces thereof, and lay 
out the vertices that will be used for the recomputation of the five-point amplitude and the 
computation of the six- and seven-point amplitudes. In section V, we display the recompu- 
tation of the five-point amplitude in great detail. In section VI, we compute and quote the 
six-point amplitude, and present the diagrams for the seven-point amplitude. We then give 
our conclusions. 

II. NOTATION 

In this section we summarize the notation used in the remainder of the paper, following 
the notation of our previous papers [11, 12]. We use the spinor helicity formalism [48, 49], 
in which the amplitudes are expressed in terms of spinor inner-products, 

(jO = = u-{kj)u+{ki) , \jl] = = u+{kj)u-{ki) , (2.1) 

where u±{k) is a massless Weyl spinor with momentum k and positive or negative chirality. 
We follow the convention that all legs are outgoing. The notation used here follows the 
standard QCD literature, with [ij] — s\^{k^k^) {ji)* so that, 

{ij)\ji]=2ki-kj = Sij. (2.2) 

These spinors are connected to Penrose's twistors [50] via a Fourier transform of half the 
variables, e.g. the U- spinors [29, 50]. (Note that the QCD- literature square bracket [ij] 
employed here differs by an overall sign compared to the notation commonly found in twistor- 
space studies [29].) We also define, as in the twistor-string literature, 

\i = u+{ki), \i = u_{ki). (2.3) 

We denote the sums of cyclicly-consecutive external momenta by 

Kt..^ = kf + kf^, + ... + k^_, + k^ , (2.4) 

where all indices are mod n for an n-gluon amplitude. The invariant mass of this vector is 
Si...j — Kl„j. Special cases include the two- and three-particle invariant masses, which are 
denoted by 

Sij = Kj^ j = (A)j -|- kj^ — 2ki ■ kj, Sij^ = (/cj -|- kj -\- kk) ■ (2-5) 
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We also define spinor strings, 

(r I (a + b) \r) = {i a) [aj] + {i b) [bj] , 

{i^\{a + b){c + d)\j-) = [ia]{a-\{c + d)\r) + [ib]{b-\{c + d)\r) , (2.6) 

and gamma matrix traces, 

ti+labcd] = [ab]{bc)[cd]{da) , (2.7) 
tr+[abc{d + e)] = [ab] {be) [cd] {da) + [ab] {be) [ce] {ea) . (2.8) 

We use the trace-based color decomposition of amplitudes [49, 51, 52, 53]. For tree-level 
amplitudes with n external gluons, this decomposition is, 

AT%{ki,hi,a,})^g^-^ ^ Tr(r"<^w • • • TM^)) ^tree(^(l'^i, . . . , ^/^n)) ^ (2.9) 

aeSn/Zn 

where g is the QCD coupling, Sn/Zn is the group of non-cyclic permutations on n symbols, 
and j^^ denotes the j-th gluon, with momentum kj, helicity hj, and adjoint color index aj. 
The T" are SU(A^c) color matrices in the fundamental representation, normalized so that 
Tr (T^T**) = 5"''. The color-ordered amplitude A^^"^^ is invariant under a cyclic permutation 
of its arguments. 

When all internal particles transform in the adjoint representation of SU(A^c), the color 
decomposition for one-loop n-gluon amplitudes is given by [54], 

K2J+1 

At^"''^\{k,,h,,a,})^g^J2''^ E E Gr„;c(<7) , (2.10) 

J C=l aeSn/Sn;c 

where [x\ is the largest integer less than or equal to x. The sum J e {0,1/2,1} runs over 
all spins of particles and nj is the multiplicity of each spin. We assume all particles are 
massless. The leading-color structure, 

Gr„;i(l) = N, Tr(r"i • • -T"") , (2.11) 

is Nc times the tree color structure. The subleading-color structures are given by 

Gr„;c(l) = Tr(r"i • --T"^-^) Tr(r"^ • • -T"") . (2.12) 

In eq. (2.10), S'„ is the set of all permutations of n objects, and Sn;c is the subset leaving 
Gr„.c invariant. For adjoint particles circulating in the loop, the subleading-color partial 
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amplitudes, An-c for c > 1, are given by a sum over permutations of the leading-cofor 
ones [3]. Therefore we need to compute directly only the leading-color, single-trace, partial 
amplitudes An-i{l^^, . . . ,71^*"). The leading-color amplitude is again invariant under cyclic 
permutation of its arguments. 

In QCD, of course, there can be fundamental representation quarks present in the loop. 
In this case only the single-trace color structure contributes, but it is smaller by a factor Nf.. 
In general, scalars or fermions in the Nc + Nc representation give a contribution, 

A'r'{{h,Ka,})^g- E ^ E Gr„;i(a)4f,(a), (2.13) 

J=0,l/2 aeSn/Z„ 

to the one-loop amplitude. We use a supersymmetric convention in which the number of 
states for a single complex scalar (squark) is 4iVc, in order to match the number of states of 
a Dirac fermion (quark). 

Helicity amplitudes that do not vanish at tree level develop infrared and ultraviolet di- 
vergences at one loop. We regulate these dimensionally. Following ref . [47] , for the divergent 
one-loop amplitudes we write, 

- cr{v^AT^ + zF^), (2.14) 
= -cr {{V^ + Kl^r + i{Ff + F^)) , (2.15) 
<i = cr ((Kf + + V:)AT + ^m + F„^)) , (2.16) 



where 



c- ^ r(i + 6)r^(i-6) 
- (47r)2- r(i - 26) ■ ^^-^^^ 

The parts contain the divergences, while the are finite. (Of course, there is some 
ambiguity in the separation between divergent and finite terms.) These pieces have a natural 
interpretation in terms of supersymmetric and non-supersymmetric parts [47], 

<r' = crA^Kf, (2.18) 
= -cr {A'rvi + iFl) , (2.19) 



cr^A^'^'V^ + zF^y (2.20) 



Here A^f ^ sums over the contributions of an A/" = 4 multiplet consisting of one gluon, 
four Majorana fermions, and three complex scalars, all in the adjoint representation. The 
J\f — 1 amplitude contains the contributions of an adjoint chiral multiplet, consisting of one 

9 



complex scalar and one Weyl fermion. The non-supersymmetric amplitudes, denoted by 



The utility of separating QCD amplitudes into supersymmetric and non-supersymmetric 
pieces follows from their differing analytic properties. The supersymmetric pieces can be 
constructed completely from four-dimensional unitarity cuts [3, 4] and have no remaining 
rational contributions (in the limit e — > 0). The examples discussed in this paper, n-gluon 
amplitudes with two negative- helicity gluons, have been known for quite some time [3, 4]. 
(Such amplitudes are often referred to as 'MHV in the supersymmetric cases, because 
amplitudes with fewer (zero or one) negative-helicity gluons vanish.) The logarithmic and 
polylogarithmic terms in the non-supersymmetric (scalar) pieces can also be obtained from 
four-dimensional cuts, or from MHV vertices. These cut-containing terms are also known 
for all n-gluon amplitudes with two negative-helicity gluons. They were computed first for 
the case where the two negative-helicity gluons are color- adjacent [4], and more recently for 
the general case [8]. 

Here we focus on the unsolved problem of computing the rational-function terms in 
the J\f — contributions in six and higher-point amplitudes, given the knowledge of the 
logarithmic and polylogarithmic terms. Finding an effective computational approach to the 
rational terms in A^.^^ is tantamount to solving the problem in QCD. 



The leading-color QCD amphtudes are expressible in terms of the different components 
in eqs. (2.18)-(2.20) via 



where nf is the number of active quark flavors in QCD. We will present the formulae for 
unrenormalized amplitudes. To carry out an MS subtraction, one should subtract from the 
leading-color partial amplitudes A^-i the quantity 



where we also included a term proportional to the number of active fundamental represen- 
tation scalars Ug, which vanishes in QCD. 

We will need to consider additional objects (parts of amplitudes), beyond the V'^ and 
defined here, in order to construct and apply appropriate on-shell recursion relations for 



Af = 0, arc just the contributions of a complex scalar in the loop, ^4^^ = A: 



{^'nv: + y/) + i(F^ + Ff)) 



(2.21) 




(2.22) 
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one-loop amplitudes. The definition of these objects and derivation of the relations is the 
subject of the next section. 

III. ON-SHELL RECURSION RELATIONS FOR LOOP AMPLITUDES 

On-shell recursion relations provide an effective means for obtaining remarkably compact 
forms for tree-level amplitudes [9, 10, 41], and have led to a variety of new results [27, 42]. 
In previous work [11, 12], we have shown how to use on-shell recursion relations to compute 
all the finite loop amplitudes of QCD. These special helicity amplitudes vanish at tree level. 
Hence the one-loop amplitudes are free of infrared and ultraviolet divergences, and they are 
'tree-like' in that they contain no cuts (absorptive parts) in four dimensions. The derivation 
of these loop recursion relations is similar in spirit to the tree-level case, but it does require 
the treatment of factorizations which differ from the 'ordinary' factorization in real momenta. 

In this paper, we will extend the analysis of refs. [11, 12] to cut-containing one- loop 
amplitudes (for which the corresponding tree- level amplitudes do not vanish), deriving new 
recursion relations for the rational functions appearing in such amplitudes. The new re- 
cursion relations allow us to systematize the factorization bootstrap approach of ref. [2]. 
We assume that the cut-containing terms have already been determined via the unitarity 
method or some other means. 

A. Analytic behavior of shifted loop amplitudes 

The starting point for our analysis, as for the finite loop amplitudes, is to consider [10] a 
complex- valued shift of the momentum of a pair of external particles in an n-point amplitude, 
kj — > kj{z), ki —>■ ki{z). This shift is best described in terms of the spinor variables A and A 
defined in eq. (2.3), 

~ ' Xi^ Xi + zXj . (3.1) 

This (j, /) shift maintains overall momentum conservation, because 

+ f^i ^ XjXj + XiXi l:j + l:i^Xj{Xj-zXi) + {Xi + zXj)Xi^^j + i(:i, (3.2) 

as well as the masslessness of the external momenta, k'j = kf = 0. Denote the original 
n-point amplitude by = An{0), and the shifted one by An{z). We wish to determine 
^n(O) by making use of the analytic properties of An{z). 
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FIG. 1: A configuration of poles and branch cuts for a term in a one- loop amplitude. The contour 
C is a circle at oo. 



In the case of tree-level or finite one- loop amplitudes, An{z) is a meromorphic function 
of z. Here we also encounter branch cuts, which may terminate at poles, as depicted in 
fig. 1. Branch cuts arise from logarithms or polylogarithms in the amplitudes. Consider, 
for example, the scalar contributions to the five-gluon amplitude with color-ordered helicity 

assignment ( recalled in eq. (4.20). It contains a logarithm, ln((— S23)/(— S51)), 

multiplied by a rational coefficient. If we perform a shift (3.1) with (j, /) = (1,2), the 
logarithm becomes 

/-..,-^(l-|3|2-h _ ,,/|231((23) + .(13m 

\-s,i + z{r\5\2-)J V(15)(|151-s|25])/ ^ ' 

This function has two branch cuts in z, one starting at 

. [1 5] 



[25]' 

the other starting at 

(2 3) 



(3.4) 



(3.5) 



(13)- 

Because of the form of the rational coefficient of ln((— S23)/(— S51)) in this case, neither 
branch cut starts at a pole. 

We assume that j and / can be chosen so that An{z) — > as z — > cx). We consider the 
following quantity, 

^{-Ar^iz), (3.6) 
2m Jc z 

where the contour integral is taken around the circle at cxd. A typical configuration for a 
term in a one-loop amplitude is shown in fig. 1. Even though the contour crosses branch 
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FIG. 2: A configuration of poles and branch cuts for a term in a one-loop amplitude, with a 
branch-cut-hugging contour. 



cuts, the integral still vanishes, because An{z) vanishes at infinity. Unlike the rational cases 
studied previously, however, this does not mean that it is given simply by a sum of residues 
at its poles. We need to include those contributions, of course; but we also need to integrate 
around the branch cuts, with special handling for poles at the end of branch cuts. 

We start with the 'ordinary' branch cuts, with no pole touching the branch cut. We can 
imagine a related contour, going along the circle at infinity, but avoiding the branch cuts by 
integrating inwards along one side, and then outwards along the other, as shown in fig. 2. 
(We will route the branch cuts so that no two overlap.) The integral along this contour 
is given by the sum of residues. The difference between the two integrals is given by the 
branch- cut-hugging integral, 

where is directed from an endpoint Bq to infinity, and B^ is directed in the opposite way. 
Now, An{z) has a branch cut along 5, which means that it has a non- vanishing discontinuity, 

27111^1^03 An{z) = An{z + ie) — An{z — ie), z on B . (3.8) 

Thus our original vanishing integral can be written as follows, 

A (z) f°° dz 

= A„,{0)+ > Res^^+ / — DiscB A„,{z) . (3.9) 

^ ' Z = Za Z J Rn Z 

poles a 

The case with a pole at the end of the branch cut — arising, for example from terms 
containing ln{—Sab)/ {^b) — can be handled similarly, but care must be taken with the 
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evaluation of the integral along the branch cut. From a conceptual point of view, we can 
also compute this case by moving the pole away from the end of the branch cut by an 
amount S, computing the branch-cut- hugging and residue terms separately, and taking the 
hmit 5 — > at the end. 



B. Cut-containing terms, and their 'completion' 

To proceed further, let us assume that we have already computed all terms having branch 
cuts, plus certain closely related terms that can generally be obtained from the same compu- 
tation. That is, we have computed all polylog terms, all log terms, and all tt^ terms. There 
are also certain classes of rational terms that are natural to include with the cut-containing 
terms. 

In particular, there are rational terms whose presence is required to cancel spurious 
singularities in the logarithmic terms. Spurious singularities arise in the course of integral 
reductions. They cannot be singularities of the final amplitude, because they are unphysical, 
and not singularities of any Feynman diagram. A simple example comes from a 'two- 
mass' triangle integral for which two of the three external legs are off-shell (massive), with 
momentum invariants si and S2, say. When there are sufficiently many loop momenta 
inserted in the numerator of this integral, it gives rise to functions such as, 

(l-r)2' ^^-^"^ 

where r is a ratio of momentum invariants (here r = S1/S2). The limit r —> 1 (that is, 
si — > S2) is a spurious singularity; it does not correspond to any physical factorization. 
Indeed, this function always shows up in the amplitude together with appropriate rational 
pieces, 

i^^^;^. (3-11) 

in a combination which is finite as r ^ 1. From a practical point of view, it is most convenient 
to 'complete' the unitarity-derived answer for the cuts by replacing functions like eq. (3.10) 
with non-singular combinations like eq. (3.11). Such completions are of course not unique; 
one could add additional rational terms free of spurious singularities. 

There are other kinds of spurious singularities connected with polylogarithms. For exam- 
ple, in the scalar contributions to the five-gluon ( — h — -|- -I-) amplitude, there are factors 
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of (2 4) and (2 5) appearing in the denominators of certain coefficients. These might appear 
to give rise to non- adjacent colfinear singularities in complex momenta; but by expanding 
the poly logarithms and logarithms in that limit, one can show that these singularities are 
in fact absent. 

Let us accordingly define two decompositions of the amplitude. The first is into 'pure-cut' 
and 'rational' pieces. The rational parts are defined by setting all logarithms, poly logarithms, 
and TT^ terms to zero, 

1 1 

(3.12) 



Rn{z) = —Ar. 

Cr 



-La 

— 

rat Cp 



ln,Li,7r2-^0 

(Note that the normalization constant cr, defined in eq. (2.17), plays no essential role in the 
following arguments, and is just carried along for completeness.) The 'pure-cut' terms are 
the remaining terms, all of which must contain logarithms, polylogarithms, or tt^ terms. 



Cn{z) = —Ar, 

Cr 



1 



Ar. 



ln,Li,7r^ 



(3.13) 



(3.14) 



pure— cut Cr 

In other words, 

Ar^{z) = Cr C„(z) + Rr^{z) 

where we have explicitly taken the ubiquitous one-loop factor cr outside of Cn{z) and Rn{z). 

The second decomposition uses the 'completed-cut' terms, obtained from Cn{z) by replac- 
ing logarithms and polylogarithms by corresponding functions free of spurious singularities. 
We shall call this completion Cn. The decomposition defines the remaining rational pieces 

Rn: 



Ar^{z) = Cr Cr^{z) + Rr^{z) 



(3.15) 



We also need to define the rational part of the completed-cut terms, CRn{z). We write. 



C„(^) = C„(z) + Ci?„(^), 



(3.16) 



where 



CRn{z) = Cr,{z) 



rat 



(3.17) 



Combining eqs. (3.14), (3.15), and (3.16), we see that the full rational part is the sum of the 
rational part of the completed-cut terms, and the remaining rational pieces. 



Rn{z) = CRn{z) + Rr,{z) . 



(3.18) 
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Now, because we know all the terms containing branch cuts, we could compute the 
branch- cut-hugging integral, 

(3.19) 



poo 
J Bo 



°° dz - 
— DiscB C„{z). 



z 



However, there is no need to do the integral explicitly, because we already know the answer 

for the integral, phis the corresponding residues. It is just C„(0), part of the final answer. 
That is, applying the same logic to Cn{z) as was applied to An(z) in eq. (3.9), we have. 



Cn{Q) = - V Res ^ - Discs C^iz) 

poles a " 



(3.20) 



In arriving at this result, we need to assume that Cn{z), in addition to the full amplitude, 
also vanishes as 2; — > 00. (This constraint may place some restrictions on the allowable 
rational completions of the cut terms.) 

Using cq. (3.9), the split-up (3.15), and cq. (3.20) to evaluate the terms involving Cn{z), 
we can write our desired answer as follows. 



-Cr 



r - Discs C^z) + y Res ^ + y Res ^ 

Jbo ^ 1 ^=^« ^ , ^=^« ^ 

■ " poles a poles a 



Cn{0)~ V] Res 



Rn{z) 



poles a 



z=Za Z 



(3.21) 



Because we have completed the cut terms so that Cn{z) contains no spurious singularities, 
the sums over the poles in eq. (3.21) are only over the genuine, 'physical' poles in the 
amplitude. (As explained elsewhere, these are the poles that arise for complex momenta, 
and not merely those that arise for real momenta.) 



C. Separate factorization of pure-cut and rational terms 

The residues of the completed-cut terms Cn{z) at the genuine poles contain both rational 
and cut-containing functions. The residues of cut-containing functions necessarily have cut- 
containing functions (and the residues of 7r^-containing terms will necessarily have factors 

of TT^, etc.), so that they will arise from cut-containing parts of the factorized amplitudes. 
The intuition from collinear factorization of one-loop amplitudes suggests, however, that 
the pure-cut terms arise from pure-cut terms, and the rational terms arise from rational 
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terms. The purpose of this subsection is to flesh out this intuition of separate factorization 
for pure-cut and rational terms. 

More concretely, the arguments of the logarithms or polylogarithms are ratios of invariants 
Si...j. (In a limited number of logarithms, the arguments are ratios of invariants to the 
renormalization scale squared.) When we shift these ratios by shifting momenta according 
to eq. (3.1), one of three things can happen to any specific argument: 

1. The ratio may be invariant under the shift, so that the residue simply has the original 
cut-containing function in it; 

2. the ratio may acquire a dependence on z, but neither vanish nor diverge at any of the 
poles in z of the given term; 

3. or the ratio may acquire a dependence on z, and either vanish or diverge at one of the 
poles in the given term. 

Because the branch cuts in a massless theory start at either a vanishing ratio or at a vanishing 
of its inverse, the last situation corresponds to a pole touching the end of a branch cut. 

In the first of these three cases, the cut-containing function clearly arises from a cut- 
containing function in the residue of the pole; and any associated rational terms arise from 
rational terms in the residue. In the second case, behavior under factorization in real mo- 
menta suggests that we need consider only single poles. In particular, for multiparticle in- 
variants, even for complex momenta we have only single poles to consider. For two-particle 
channels, double poles can arise (for complex momenta) in the kinematic invariant of the 
momenta of two nearest-neighboring legs in the color ordering, but only at one loop and 
only for certain helicity configurations. 

In two-particle factorizations, one-loop splitting amplitudes appear in addition to tree- 
level ones. The remaining amplitude left behind in the former case, however, is a tree-level 
amplitude, which is of course purely rational. In the latter case, we need consider only single 
poles. 

We are left to consider the separation inside the cut-containing one-loop splitting ampli- 
tudes themselves. Because we are only interested in the scalar loop contributions, the rele- 
vant sphtting amphtudes we need are scalar-loop ones [3, 4, 15, 28, 55]. For the helicity con- 
figurations for which the tree- level splitting amplitude Split ^^'^ (a; ; a'*'", 6'*''') is non- vanishing 
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(where A is the hehcity of the outgoing off-shell leg, and x is the longitudinal-momentum 
fraction carried by leg a), we write the loop splitting amplitude as, 

gp^.^i-ioop, [0] ^A. ^ ^A,) _ ^1-ioop A„A„ [0] ^^^^ ^ Split!!^^(a;; a^" , 6^") , (3.22) 

where 



^5 



i-ioop-+,[o]^^^^^ = 0, (3.23) 

i-ioop++,[o]^^ _ (l^^Y 2.t(1-3;) 

s) (l-2e)(2-2e)(3-2e) 



■s {x, 5) = cr ( ^ ) . (3.24) 



The vertex corresponding to the first splitting amplitude will indeed vanish, and so it will 
not affect the separation between cut-containing and rational pieces. The situation with 
the second helicity configuration remains to be studied. Indeed, it appears in this case that 
for complex momenta, there are genuinely non-factorizing contributions in addition to pole 
contributions. The shift choices we will make in this paper avoid the appearance of this 
vertex, and so we postpone the analysis of this case to future work. 

Of course, outside of cut-containing terms, double poles do arise in general, as indicated 
by the finite one-loop splitting amplitude (the helicity configuration that vanishes at tree 
level) , 

Splitr°°^'™(a;;a+,6+) = -^.V<^)^ ■ (3-25) 

They would be handled as in refs. [11, 12]; given the shift choices we will make in the present 
paper, they will in any event not arise. 

So long as we are indeed considering only single poles, extracting the residue will not force 
us to expand C„(z) in a series around the pole. In this case, the two types of contributions, 
pure-cut and rational, will remain separate in this class of contributions. As explained above, 
the third and last case, that of the pole hitting the branch cut, can effectively be reduced 
to the second case by artificially separating the pole from the end of the branch cut, and 
removing the separation at the end of the calculation. 

As an example of the different kinds of behavior, consider the following expression. 



(3.26) 



(2 3) [34]^ (41)^ (12) [15]Li(^rfg 
(3 4) (4 5) si 



where Li(r) = (ln(r) + 1 — r)/(l — r)^, under a (1, 5) and (3, 4) shift respectively. Under the 
first shift, we obtain a simple pole at ^ = — (4 5) / (4 1), but the argument of the Li function 
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is unchanged, and accordingly, the logarithm and rational parts are left unaltered. In the 
context of a factorization, wc would interpret the logarithm as arising from a logarithm in 
the factorized amplitude, and the rational part inside the Li from a rational term. 

Under the second shift, we obtain a simple pole at 2; = — (4 5) / (3 5). The argument of 
Li at the pole is now, 

-(2 3)([3 2]-;.poie[4 2]) _ (2 3) ((5 3) [3 2] + (5 4) [4 2]) _ (2 3) (1 5) [1 2] ^ 
-S51 -(3 5)s5i -(3 5)s5i 

but because the pole is a simple one, the separation between logarithm and rational part 
is again left undisturbed. If instead the original expression had contained not a (4 5) in 
the denominator but rather a (4 5)^, we would have had to expand the logarithm inside Li 
in order to extract the residue, and this expansion would have led to rational terms whose 
origin was in a 'pure-cut' term, rather than a rational one, mixing up the two types of 
contributions. Fortunately, such combinations never arise in the analysis of the amplitudes 
discussed in this paper. 



D. Residues of the remaining rational pieces Rn{z) 

Thanks to the analysis of the previous subsection, we can simplify our calculation by 
separating the two classes of terms — pure-cut and rational — in the factorized amplitudes. 
Because we already know the cut-containing pieces, we need to analyze only the rational 
terms. (This separation also avoids the need to treat poles hitting the ends of branch cuts 
explicitly.) 

We now examine more carefully the residues of the poles of the purely rational terms, 
Rn{z), and related to that, of the shifted remaining rational terms, Rn{z), appearing in 
eq. (3.21). Given the (j, Z) shift (3.1), we define a partition P to be a set of two or more 
cyclicly-consecutive momentum labels containing j, such that the complementary set P 
consists of two or more cyclicly-consecutive labels containing /: 

P = {Pi,P2,...,i,...,P-i}, (3.28) 
P = {Pi, P2,..., P-i}, 
P U P = {l,2,...,n}. 

This definition ensures that the sum of momenta in each partition is ^-dependent, so that 
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it can go on shell for a suitable value of z. At tree level, the sum over residues becomes, 



Z — Zoc 

poles a 



partitions P h=± 

X p2 '2 X ^Tikp,, ... i/, ... , P-'^) . (3.29) 

1 III p 



The complex on-shell momenta kj, ki and P are determined by solving the on-shell condition, 
P^ = mp, for 2;. Although the examples we discuss in this paper are for massless particles, 
we allow for a mass mp to indicate that there is no special restriction to massless amplitudes, 
as has already been noted at tree level [27]. 

At one loop, the sum analogous to eq. (3.29) will have an additional two-fold sum. In 
each term in this sum, either Al or Ar will be a tree amplitude, and the other one will 
be a loop amplitude; in general both terms will appear. Taking the rational parts of the 
one-loop amplitudes appearing in this expression, the one-loop physical-pole recursion for 
the rational terms is, 

- V Res^ ^ R^{k,,...,kn) 



Z=Za Z 

poles a 



partitions P /i=± 



X p2 X ^ (^Pi I • • • I ^i) • • • ) ^p_i ) P ) 

+ ^*'-^^(A;p„..., A;p_„-P'^) (3.30) 

X R{kp^ , . . . ,ki, . . . , kp_^ , P '*) , 



^P^ 



where we now assume that the intermediate states are massless. This result follows directly 
from the general factorization behavior of one-loop amplitudes, plus the separate factoriza- 
tion of pure-cut and rational terms that was established in the previous subsection. Just as 
in the case of the tree-level recursion (3.29), it exhibits the required factorization properties 
in each channel P (dropping the terms with logarithms, polylogarithms, and tt^). Although 
the R functions are not complete amplitudes, they can be thought of as vertices from a 
diagrammatic perspective, and this equation lends itself to the same kind of diagrammatic 
interpretation available for eq. (3.29). 
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However, the factorization cannot distinguish between the rational terms we have already 
included in the completed-cut terms and the remaining ones. That is, there would be an 
overlap or double count if we were simply to combine the recursive diagrams with the 
completed-cut terms. To remove this overlap, we separate the physical-pole contributions 
into those already included in the completed-cut terms and those in the remaining rational 
terms. Using eq. (3.18), we know that 

- y Res^ = i?^ = - y Res^^- V Res ^ . (3.31) 

^ Z=Zc, Z ^ Z=Zc Z ^ Z=Zoc Z 

poles a. poles a poles a 

Because we know the completed-cut terms Cn{z) and their rational parts CRn{z) explicitly, 
we can compute the first term on the right-hand side, and solve for the remaining terms 
using the determination of via eq. (3.30). Inserting the result into eq. (3.21) then gives 
us the basic on-shell recursion relation for complete one-loop amplitudes, 

An(0) = cr \Cn(0) + R^+ y Res ^iM^] . (3.32) 

' Z=Za Z J 

poles a 

To compute with this equation, we construct R^ via 'direct recursion' diagrams; that is, 
via eq. (3.30). We call the elements of the last term 'overlap' terms. Because each pole is 
associated with a specific diagram, these can also be given a diagrammatic interpretation. 
Although the definition of the completed-cut terms C„ is not unique, the ambiguity cancels 
between Cn(0) and the sum over Ci?„ residues. In the calculations in this paper, an astute 
choice of completion terms can simplify the calculation, by simplifying the extraction of the 
residues in the last term. 

The reader may wonder how the calculation would have proceeded if we had started with 
'pure' cut terms, not including any of the rational pieces needed to eliminate the spurious 
singularities. In this case, the intermediate stages, and in particular eq. (3.21), would have 
included a sum over the spurious singularities as well. Since these singularities include 
double and triple poles, we would have needed to expand the logarithms in extracting the 
residues for the 'overlap' terms, and this expansion would have produced rational terms. 
In contrast, with our approach, we never evaluate residues at values of z corresponding to 
unphysical spurious singularities, and we never have to expand logarithmic functions. 

In section II, we separated the one-loop amplitudes into divergent and finite parts. The 
amplitude as a whole satisfies the bootstrap relation eq. (3.32); but it turns out that for the 
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amplitudes we consider in the present paper, it can be applied separately to the V and F 
terms. As the recursion relation for the pure-scalar parts of the former are basically the same 
as at tree level, we will focus on the computation of the F terms. For this purpose, we shall 
use quantities analogous to C„, Cn, Rn, and CRn, defined as in eqs. (3.13)-(3.17), but 
with respect to F^ of eq. (2.14) instead of An- Note that this shift of convention generates 
a relative factor of i in the quantities we use below, due to the relative i in eq. (2.14). 

IV. REVIEW OF KNOWN RESULTS 

In this section we summarize the previously- computed amplitudes, and pieces thereof, 
that feed into our recursive construction. In this paper, we consider n-gluon amplitudes with 
two color- adjacent negative helicities, A„.i(l~, 2~, 3+, 4+, . . . , (n — 1)+, n'^). As mentioned 
in section II, the — A and M —1 components of these amphtudes have been known for a 
while, so the only issue is the computation of the jV = or scalar contribution, 

A[:!l(l-,2-,3+,4+,...,(n-l)+,n+). (4.1) 

Assigning intermediate helicities to all possible factorizations of this amplitude, as encoun- 
tered in eq. (3.30), allows us to determine which lower-point amplitudes are required as input. 
Besides the Parke- Taylor tree amplitudes with two (adjacent) negative hehcities [52, 53, 56], 
we shall need the one- loop scalar contributions with one negative helicity [57], for which 
we recently found a compact form [12]. The one-loop amplitudes with two adjacent nega- 
tive helicities and smaller values of n are also needed, and are obtained recursively, given a 
suitable starting point. 

By choosing to shift the two negative-hehcity legs we can avoid some factorizations. For 
example, for a generic choice of j and I in the (j, I) shift (3.1), the amplitude could factorize 
onto products containing a one-loop amplitude with all positive external helicities and an 
internal helicity of either sign. We avoid these factorizations by choosing to shift the two 
negative-helicity legs, (j, /) = (1,2). One external negative helicity then appears in each 
partition. In addition to lower-point amplitudes, we also need the logarithmic parts of the 
n-point amphtude (4.1), obtained in ref. [4]. (The logarithmic parts of the more general set 
of amplitudes for two non-adjacent negative helicities were recently obtained in ref. [8].) 

In section V we recompute the five-gluon helicity amplitude >l5°!^(l~, 2~, 3+, 4+, 5+). For 
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that computation, we take the four-point amphtude ^^^'^^(l^, 2^, 3^, 4"'") as an input, but 
here we quote the previous five-gluon result too, in order to demonstrate that the method 
works. (We also outline the recursive construction of A|°!l(1~, 2"*", 3~, 4+, 5"*"), which uses 
^4!i(l~> 2+, 3~, 4+) as an input, and agrees with the known result.) In section VI we feed 
the five-point result back into the recursion to construct the six-point result, and following 
that, outline the construction of the seven-point result. 

Finally, we also need the three-point amplitudes, which vanish for real momenta, but are 
non-vanishing for generic complex momenta. The one-loop three-vertex that we need may 
be deduced from a one-loop splitting amplitude [3]. 

Let us start with the tree amplitudes. The tree amplitudes that enter into our calculation 
are just the Parke- Taylor amplitudes [52, 53, 56], 

/ir(1^2+,3+,...,n+) = 0, (4.2) 

For some complex momenta, the three-point amplitudes 

(12^^ 

are non- vanishing. (They are vanishing for real momenta.) 

The finite one-loop amplitudes that feed into our recursion are also relatively simple. The 
four-point finite amplitude with a single negative-helicity leg is [58, 59, 60], 

A,.^,[l ,2 ,3 ,4 )- g [12] (2 3) (34) [41]- ^^'^^ 

Since the amplitude is entirely composed of rational functions, the vertex is proportional to 
the amplitude 

R,{r, 2+, 3+, 4+) = — 4°1(1", 2+, 3+, 4+) . (4.7) 

We have removed an extra i from this vertex (and all others in the section) compared to the 
vertices Rn appearing in section III. As mentioned at the end of that section, we wish to 
perform the recursion directly on the finite parts defined in eq. (2.14). For this reason, 
a factor of i from the vertex is removed from the vertex, compared with the one that would 
be used for constructing the amplitudes An. 
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The five-point finite amplitudes are also rather simple. We will need the finite amplitude 



4°l(l-,2+,3+,4+,5+) = 1% ^ 



[2 5]= 



+ 



(14)^ [4 5] (3 5) (13)^ [3 2] (4 2) 



3 (34)' L [12] [51] (12) (2 3) (4 5)' (1 5) (5 4) (3 2)'. 

This amplitude, along with all the other one-loop five-gluon helicity amplitudes, was first 
calculated using string-based methods [47] . Again because the amplitude is purely rational, 
the vertex is proportional to the amplitude, 

1 



i?5(l-, 2+, 3+, 4+, 5+) = — 2+, 3+, 4+, 5^ 



(4.8) 



It is worth noting that compact expressions now exist for the n-point generalization of 
this amplitude [11, 12], 74|^\(1~, 2+, 3+, . . . , n"*"), which agree numerically with Mahlon's [57] 
original determination. 

Using the decompositions (2.14)-(2.21), we express the divergent four-point amplitudes 
^•!|(l-,2-,3+,4+) in terms of the functions [58, 59, 60], 

,2 ^ f , ,,2 



V 



V-S12/ V-S23^ 

-(— 



3 



/ 



' 3 9 

In this case the finite parts are trivial. 



(4.9) 
(4.10) 
(4.11) 



= . 



(4.12) 



The version of dimensional regularization under consideration is determined by the 8r pa- 
rameter; for the 't Hooft-Veltman scheme [26] we take 5r — 1 while for the four-dimensional 
hehcity scheme [58, 61] we take 5r — 0. 

Equation (4.12) shows that all the rational terms for the ( 1-+) case arc constant 

multiples of the tree amplitude, so they can easily be absorbed into the terms. Because 
the tree amplitude obeys its own on-shell recursion relation [9], it does not really matter 
whether we put constant terms hke the 2/9 term in VI into the y or F category, but we 
will be able to drop one recursive diagram by assigning it to V . In general, then, we define 

by, 

2 

K^ = -^ + 9- (4.13) 
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Because Fl vanishes, we take the loop vertex to also vanish, 

i?4(l-,2-,3+,4+) = 0. (4.14) 

For the divergent amplitudes with five or more legs, it is useful to introduce a set of 
auxiliary functions [47], 

ln(r) 



Lo(r) 



1 — r ' 

ln(r) + 1 — r 



(1-r) 



2 



,^^,,.m^±zim, (,,5, 

in which the pole at r = 1 is removable. As discussed in section III, we can therefore 
use the functions to construct the completed-cut terms, out of logarithms deduced from 
four- dimensional cuts. 

We shall be quoting the functional form of the cut-containing pieces in the Euclidean 
region; a discussion of analytic continuations to the physical region may be found in, for 
example, ref. [2]. 

We will also need the five-gluon amplitude 745"[|(1~, 2~, 3+, 4+, 5+). Ref. [47] gives us, 
Vi = Ati^y + h<^^) ln(^^) 4-^-^, (4.16) 

2 X e , ,,2 

2, (4.17) 



^ 2e 



= + (4.18) 

for the functions appearing in the decompositions (2.14)-(2.21). (The 1/e singularities differ 
from those in ref. [47] because the amplitudes there were renormalized, whereas here we 
are using unrenormalized amplitudes; the difference is given simply by the renormalization 
subtraction (2.22).) The finite parts of this amplitude are, 

^ _ 1 (1 2)^ ((2 3) [3 4] (4 1) + (2 4) [4 5] (5 1)) Lp (efg) 

~ 2 (2 3) (3 4) (4 5) (51) ^^'^^^ 



1-/ 1 



[3 4] (4 1) (2 4) [4 5] (2 3) [3 4] (4 1) + (2 4) [4 5] (5 1) U 



-S23 



3^' 3 (34) (4 5) sii 

+ ^5 , (4.20) 
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where 



(3 5) [3 5f 



+ 



(12) [3 5f 



3 [12] [2 3] (3 4) (4 5) [51] 3 [2 3] (3 4) (4 5) [5 1] 
1 (1 2) [3 4] (41) (2 4) [4 5] 
6 S23 (3 4)(4 5)s5i 



(4.21) 



In section V, we shall recompute the explicit rational terms R^. The other pieces are all 
either trivial or obtainable from four-dimensional unitarity. 

Following the discussion of the previous section, we define a recursion vertex (3.12), 
composed of all the rational terms in F^, 



it:5(l-,2-,3+,4+,5+)-F5 



rat 



(4.22) 



including those contained in L2. As in the four-point case the effect of the 2/9 term in 
is trivial, so we do not need to make it part of the recursion vertex. 

Now consider the known six-point results for A^"[|(l~, 2~, 3+, 4+, 5+, 6+). Except for ratio- 
nal terms associated with the J = scalar loop, these amplitudes were determined from the 
unitarity method in rcfs. [3, 4], where the calculations were performed for the more general 
n-point amphtudes with two adjacent negative heli cities. These results have recently been 
confirmed in refs. [8, 22]. 

Prom the results of ref. [3], after setting n = 6, we have the J\f — A contribution, 

/_« N ^^^^^^^^2^ ^423) 



where all indices are to be taken mod 6 and 



-VLi2(l- 



1 ^ 



j^^l \ ""Si+l,i+2,i+3 

From ref. [4] , the J\f — 1 components are. 



In 



(4.24) 



-1// 2\s / 2 

Vq — 



2e V V -S23 



-sqi 



2, 



{i2y 



2su (2 3) (34) (4 5) (5 6) (61) 



- Lo (— ) (tr+[1256] - tr+[12(l + 6)5] 

'16 ^ S16 / V 

+ — Lo (—) ftr+[1234] - tr+[124(2 + 3)] 



(4.25) 
(4.26) 
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while the scalar loop contributions are, 



6 Q 




1 (12)^ 



35?2 (2 3) (3 4) (4 5) (5 6) (61) 



X ^ L2 (— ) ((tr+[1256])2 tr+[12(l + 6)5] - tr+[1256] (tr+[12(l + 6)5])^ 




+ ^L2(— ) (tr+[124(2 + 3)](tr+[1234])2 - (tr+[124(2 + 3)])^ tr+[1234] 



*234 ^'5234^ ^ 




where are the rational terms not contained in L2. A key task of this paper will be to 
obtain an explicit formula for the unknown rational function Rq. We shall do so in section VI. 

Similarly, for the seven-point amplitude ^^^^(l", 2", 3+, 4+, 5+, 6+, 7+), it is not difficult 
to extract the functions from the n-point forms given in refs. [3, 4]. Although we shall not 
discuss the seven-point case in any detail, these functions enter into the computation of the 
rational terms of this amplitude, as outlined in section VI. 

We still need a one-loop three-vertex for the recursion. We determine this vertex by 
inspecting the one-loop splitting amplitudes. Because the loop splitting amplitude with 
opposite on-shell helicities and a scalar circulating in the loop vanishes (see eq. (3.23)), 



For the cases where the two external lines have the same helicity, the situation is much more 
subtle, with the appearance of 'unreal poles' [11, 12] and non-factorizing contributions. 
However, because we choose to shift the two negative-helicity legs, we do not encounter such 
vertices. 

Finally, we need the rational functions Ci?„ contained in the completed-cut part defined 

in eqs. (3.15) and (3.17), for n = 5, 6, and 7. These functions will be used to obtain the 
overlap contributions in sections V and VI. We easily obtain these from the L2 terms by 
replacing the L2 function with its rational part, using eq. (4.15). For the five-point amplitude 



SplitL"Ul^,2^) = 0, 



(4.29) 



the corresponding three-vertex should also be taken to vanish. 



i?3(i-,2+,-^±) = o. 



(4.30) 
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/lI['.\(l",2~,3+,4+,5+), the explicit value is, 



Ci?5 — — - 



Sl5 + S23 



[34] (41) (24) [4 5] (2 3) [34] (41) + (24) [45] (5 1) 



(4.31) 



6S23Sl5(Sl5-S23)' (3 4) (4 5) 

For the six-point amplitude A|!!^(1", 2~, 3+, 4+, 5+, 6+), the rational function appearing in 
the completed-cut part is, 



CRr — — - 



1 



6 (2 3) (3 4) (4 5) (5 6) (61) 

1 SlQ + S234 



X 



(2 5) [5 6] (61) (5 1) (2-| (3 + 4) |5-> 



SI6S234 (Sl6 - -5234)^ 

X I (2 5) [5 6] (6 1) + (2- 1 (3 + 4) |5-) (5 1)) 



+ 



1 



S234 + S23 



'S23"S234 (S234 — S23J 



(2 4) (4+ 1 (5 + 6) 1 1+) (2 3) [3 4] (4 1) 



X ((2 3) [34](41) + (2 4)(4+|(5 + 6)|l+>) 



(4.32) 



Similarly, for the seven-point amplitude 747°!l(1~, 2~, 3"*", 4+, 5+, 6"*", 7"*") which we briefly dis- 
cuss in section VI, it is not difficult to extract the explicit value of CR7 from eq. (7.1) of 
ref . [4] , although we shall not quote the result here. 



V. RECOMPUTATION OF FIVE-GLUON QCD AMPLITUDES 

In this section we illustrate our method for determining loop amplitudes, by recomputing 
the known five-gluon QCD amplitudes, given the four-dimensional cut-constructible parts 
of the amplitudes. There are two independent helicity amplitudes, 74|!!^(1~, 2~, 3+, 4+, 5+) 
and Aji(l-,2+,3-,4+,5+). We will discuss the first of these in some detail, and merely 
summarize the calculation of the latter. In both cases, we correctly reproduce the results of 
ref. [47]. 

Begin with 745°[(1~, 2~, 3"*", 4+, 5+). For this amplitude, use a (1,2) shift, 

Ai — > Ai , 

Ai — > Ai — ^A2 , 

A2 — A2 + zXi , 

A2 ^ A2. (5.1) 
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3\ ^4+ 3\ 










(d) (e) ^i- ^/ (f) 

FIG. 3: The recursive diagrams for computing the rational parts of A5!'^(l~,2~,3+,4+,5+) with 
the shift of legs 1,2 given in eq. (5.1). T' signifies a tree vertex and 'L' a loop vertex. 

It is not difficult to verify that this shift has the required property that the rational part of 
the cut term CR^, given in eq. (4.31), vanishes at large as required. 

This shift yields a version of the rational-recursion (3.30), where each term is represented 
by one of the recursive diagrams depicted in fig. 3. We have dropped diagrams with a 
trivially vanishing tree amplitude. Consider the first diagram in fig. 3, 



L'f = >lf ^(2-, 3+, -X-) X — X i?4(i-, kt^. 4+, 5+) . (5.2) 

523 

It vanishes, 

Df^Q, (5.3) 

because [9] 



^r(2",3+,-if2-3) CX (2 if 23 

:2-|(2 + 3)l2-)- /_:' (1-1(2 + 3) |2 



oc 



(1-|(2 + 3)|2- 

= 0. (5.4) 
Diagram (b) also vanishes, 

Df^ = , (5.5) 

because the loop three-vertex (4.30) vanishes. Similarly, it is not difficult to show that 
diagrams (d) and (e) vanish, 

= ) = . (5.6) 
We are left with just two direct-recursion diagrams. Diagram (c) is given by 



L>f) =^*r(2-,3+,-i^2+) X — xi?4(i-,i^23,4+,5+). (5.7) 
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FIG. 4: The five-point overlap diagrams using the (1,2) shift given in eq. (5.1). 



As we saw in eq. (4.14), the loop vertex vanishes, and so 

= . (5.8) 

The last diagram is, 

L>f = Ar^(5+,i-,-K5l) X — X i?4(2-,3+,4+,K5+) 

1 {H-K5i)f 1 (3i^5l) [3i^5l]' 



3 (5i) ((-i^5i)5) 551 [2 3] (3 4) (4^15) [K,i2] 

1 (3-|4|2-) (1-|5|3-) 



1 (1- 


5 


2- 


r 


3 (51) (5" 


1 


2-) 



-\3 



(51)[15](1-|5|2-)^ 

[2 4] [3 5]^ 



(5.9) 



3(34) [12] [15] [2 3]'' 
This direct-recursion diagram is the only one that does not vanish. 

Next we must evaluate the overlap contributions from eq. (3.32), depicted in fig. 4. We 
start from the rational parts of the cut contributions, CR5 as given in eq. (4.31). Applying 
the shift (5.1) to these contributions, we have, 
1[3 4] (41)((2 4) + z(14))[4 5] 



6 (3 4) (4 5) 

X f((23) + z(13))[34](41) + ((24) + z(14))[45](51)) (5.10) 



S51 + S23 - z{l 


5 


2-) + z{l- 


3 


2-) 


((2 3) + ^ (1 3)) [3 2] (1 5) ([5 l]-z[5 2])(si5 - ^23 - 


■z{l- 


(5 + 3)|2-))2 



The residues of CR^{z)/z that we need to evaluate are located at the values of z, 

^ " (13)' [2 5]' ^^-^^^ 

corresponding to the two overlap diagrams in fig. 4. Evaluating the residue corresponding 

to the first of these overlap diagrams is very simple and gives, 

^(a) ^ 1 (12)^ (14) [3 4] 
' 6 (15) (2 3) (3 4) (4 5) [2 3] ' ^ ' ' 
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Similarly, the overlap diagram (b) gives 

Q(b) ^ 1 (14) [3 4] [3 5] ((14) [3 4] -(15) [3 5]) 
' 6 (15) (3 4) (4 5) [15] [2 3]' ' ^ ' ' 

Summing over the non- vanishing diagrammatic contributions, we get the simple result, 

^5 = + of + of) 

If 2 [2 4] [3 5]^ (12)^ (14) [3 4] 



6V (3 4) [12] [15] [2 3]' (1 5) (2 3) (3 4) (4 5) [2 3] 

(14)^ [3 4]^ [3 5] _ (14) [3 4] [3 5]^ \ 
(15) (34) (4 5) [15] [2 3]' (3 4) (4 5) [1 5] [2 3] V ' ^' ' 

With a few spinor manipulations this result can be brought into manifest agreement 
with the known result (4.21). Thus, we have correctly reproduced the rational parts of 
74|!!^(1~, 2~, 3+, 4+, 5+), without performing any loop integrals. 

We have also verified that our method properly reproduces ^^^^^^(l", 2+, 3~, 4+, 5+), com- 
puted in rcf. [47]. As for the previous case, we choose the negative- helicity legs as the shifted 
ones, i.e. we use a (1,3) shift, 

Ai — > Ai , 

Ai ^ Ai — ZX3 , 

A3 — > A3 -|- zXi , 

A3 ^ A3. (5.15) 

For this computation, we completed the cut terms with the Lsi function of ref. [47], along 
with the Li and L2 functions also used there, 

(12)(23)(34)(41)'[24]' 2 Ls,(^, ^) + L, (^) + L, (^) 



(4 5) (51) (2 4)' si 



(32)(21)(15)(53) 
(5 4) (4 3) (2 5)' 



[2 5]'2Ls,(5f^,5f^)+L,(5f^)+L,(5|t) 



'34 



2 (23)'(41)^[24]^L,(^eftJ 2 (2 1)^(5 3)^ [2 5]^ (^^f^^ 
^ 3 (45)(51)(24) si 3 (54)(43)(25) ^ " ^ 

L2 (^) .1 (13) [2 4] [2 5] ((1 5) [5 2] (2 3) - (3 4) [4 2] (2 1)) 



+ si V3 (4 5) 



2(12)^(3 4)' (41) [2 4]-'^ _ 2(3 2)^(15)^(5 3) [2 5]^ 
^ 3 (4 5) (5 1) (2 4) 3 (5 4) (4 3) (2 5) 
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(a), (b), (d) and (e) vanish. 

1 (1 3)=^((1 5) [5 2] (2 3) - (3 4) [4 2] (2 1)) Lq (fft) 
6 (12) (2 3) (3 4) (4 5) (51) S51 

This function satisfies C^iz) — > as z ^ 00, thanks to cancellations between the polyloga- 
rithmic or logarithmic functions and the rational terms. The attentive reader will note that 
this function has spurious double poles involving non-adjacent legs in the color ordering, 
e.g. 1/ (2 4)^. These poles do not invahdate the calculation with a (1,3) shift, because 
they acquire no z dependence, and hence produce no poles in z at spurious locations. The 
spurious singularities cancel in the complete answer for F^, not only at order 1/ (2 4)^, but 
also at order 1/ (2 4), even for complex momenta. (The completed-cut term C5 given here 
does contain a 1/ (2 4) pole for complex momenta, but it is cancelled in the full answer by 
the additional diagrammatic terms in our construction.) 

VI. SIX- AND SEVEN-POINT QCD AMPLITUDES 

In this section we describe the computations of the unknown rational func- 
tions for the six- and seven-point helicity amplitudes Ag^!L(l~, 2", 3+, 4+, 5+, 6+) and 
74y°[(l~, 2~, 3"*", 4+, 5"*", 6"*", 7"*"). Were one to attempt the calculation by traditional means, 
one would encounter large numbers of Feynman diagrams. The total number of one-loop 
diagrams for the six-gluon process in QCD is 10,860 (including gluon, ghost and fermion 
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loops, but dropping those which vanish trivially in dimensional regularization) . For seven 
and eight external gluons the numbers grow to 168,925 and 3,017,490 respectively. But the 
number of diagrams only hints at the full complexity of the calculation, because it does not 
take into account the explosion of terms resulting from tensor integral reductions, which is 
what renders a brute-force Feynman diagram computation impractical even on a modern 
computer. The approach we take in the present paper avoids this explosion by focusing on 
analytic properties that all amplitudes must satisfy. We computed the logarithmic terms 
in the amplitudes long ago [3, 4], using the unitarity-based method. Here we complete 
the QCD calculation by computing the rational terms in the scalar contributions, namely 
eq. (4.1) for n — 6,7. For the six-point case we present a compact analytical expression. 

First consider the six-point amplitude 74g°!^(l~, 2', 3+, 4+, 5+, 6+). To obtain the rational 
terms of the six-point amplitude, Rq in eq. (4.28), we first evaluate the recursive diagrams 
shown in fig. 5, corresponding to the terms in the recursion (3.30). We again shift the two 
negative-helicity legs, using cq. (5.1). Following a similar discussion as for the five-point 
amplitude, it is not difficult to show that diagrams (a), (b), (d), and (e) vanish, 

= D^^^ = L>f = D^^'^ = . (6.1) 

Four diagrams remain to be evaluated, 

D, = L>(^) + + L>f + Di"^ 

= Ar(2-,3+,-i^+3) X — X i?5(i-,i^23,4+,5+,6+) 

+ R,i2-, 3+, 4+, 5+, K+) X — X A*r(6+, i", -K^,) 

+ R,{2-, 3+, 4+, -^+4) X — X ^*-^(i-, ^234, 5+, 6+) 

S234 

+ A'n-K,,„ 2-, 3+, 4+) X — X i?4(i-, Kts4, 5+, 6+) . (6.2) 

•§234 

We do not present a detailed evaluation of these diagrams, because it is similar to the 
five-point evaluation discussed in the previous section. 

Next consider the overlap contributions of eq. (3.32) displayed in fig. 6. These are de- 
termined by evaluating the residues of the rational part of the cut contributions given in 
eq. (4.32) after performing the shift (5.1). The residues of CRq{z)/ z that need to be com- 
puted are at the values of z located at, 

[2 6]' (13)' (1-| (3 + 4) |2-) ' ^ ' 
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FIG. 6: The overlap diagrams for 2", 3+, 4+, 5+, 6+), using the (1,2) shift. 

corresponding to the three overlap diagrams in fig. 6. It is straightforward to evaluate the 
residues since we encounter only simple poles. For diagram (a) of fig. 6, for example, we 

have, 

Q(a) ^ 1 (1 5) [5 6] (5-| (1 + 2) |6-) (2(15) ^612 - (11 2(1 + 6) |5+)) 
6 (34)(45)(56)(61)[61](3-|(l + 2)|6-)(5-|(3 + 4)|2-)^ ' 

which has been simplified using spinor-product identities. Similarly, it is not difficult to 
evaluate the other two residues, corresponding to diagrams (b) and (c) in fig. 6. 

The sum of the non- vanishing recursive and overlap diagrams for in eq. (4.28) is given 

by 

^6 = + + ) + dS"^ + Ot^ + Of^ + 0^^ . (6.5) 

One of the terms in the expression for the overlap diagram in the (2 3) channel, 0^^\ contains 
an unusual factor in its denominator, namely the square of the quantity 

(l-|5(l + 2) |3+) + (13)s56 = (r| (2 + 3)4|3+> - (r|6(l + 2) |3+> . (6.6) 

However, the recursive diagram in the (2 3) channel, D^^\ shown in fig. 5(c), also contains a 
term with this behavior, and the two terms cancel against each other. After simplifying the 

sum over diagrams we obtain, 

R ^ lf_2(35)[35](4-|(l + 2)|3-)(4-|(l + 2)|6-)(5-|(l + 2)|6-) 



6 1 [12] (3 4)14 5)' [6 1] (5"! (3 + 4) 12") (6"! (1 + 2) 13" 
(3 5) [3 6] (4-|(l + 2)|6-)^ 
[12] (34)^(45)161] (51(3 + 4)12-) 

(12) (24) (35) [35]l56] (51 (1 + 2) |61 



- 2: 

+ 2 
+ 2 



(34)' (45) [6 1] (21 (1 + 6) |5-) (51 (3 + 4) |21 (61 (1 + 2) 13") 

(1 2)^[3 5]^((5-| (3 + 4) 2 |1+) + (5-| 35 |1+)) 

(34) (45) (6 1) (21 (1 + 6) |51 (51 (3 + 4) 12") (61 (1 + 2) 13") 

(12)-^ (3 5) [4 6] [5 6] ,^ [3 6]' 



(2 3) (3 4) (4 5) (5 6) (1" | (2 + 3) 14") (3" | (1 + 2) 16") [1 2] [2 3] (4 5)^ [6 1] 
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[56] (5-| (1 + 2) |6-)'(2 (4-| (3 + 5)(1 + 2) |5+) + (1 2) [1 2] (45)) 
[1 2] (3 4) (4 5)' (5 6) [6 1] (3" | (1 + 2) |6-) (5" | (3 + 4) |2-) 
^ , (1 5)^ [3 Af [5 6] ((1 6) [3 4] (4 5) - (l" | (2 + 4) |3') (5 6)) 
[2 3] (4 5) (5 6)^^234 (1-| (2 + 3) 14") (5-| (3 + 4) 12") 



(12) (15) [3 4] [5 6 


(1- 


(5 + 6)(3 + 4) |5+) 


(3 4) (4 5) (5 6)s234 (l- 


(2 + 3)|4-)(5- 


(3 + 4)|2-) 



(3 5) (1- 


(2 + 4) |3-)' 


[2 3] (34) (4 5) (5 6) (61)s234 (5" 


(3 + 4)|2-) 



+ 2 

(1 2) (1-| (2 + 4) |3-) (2 (1-| (2 + 4) |3-) + (I'l 4 |3-)) 
[2 3] (3 4) (4 5) (5 6) (6 1) S234 

+ 2- 



+ 2 



(12)'[4 6]'(5- 


(4 + 6) 


5-) 


(2 3) (4 5) (5 6) S123 (1-|(2 + 3)|4-) (3" 
(12)^[3 5]^(4-| (3 + 5) |4-) 


(1 + 2)|6-) 



(34) (45) (6 1) S612 (2-| (1 + 6) 15") (6"! (1 + 2) 13") 



(12)' 



(2 3) (34) (4 5) (5 6) (61) 



;i-|4|3-) ^ (2-|5|6-) 



(6.7) 



[2 3] [61] 

determining the previously unknown rational function in eq. (4.28). 

As in the past, for both theoretical and practical reasons, it probably will be important 
to have the simplest representations of amplitudes. It is likely that the result in eq. (6.7) can 
be simplified even further. For example, individual terms contain spurious singularities due 
to the following denominator factors: (1"! (2 + 3) 14"), (2"! (1 + 6) 15"), (3"! (1 + 2) 16"), 
(5~| (3 + 4) |2~), and (6~| (1 + 2) |3~). We have checked numerically that these singularities 
cancel between different terms, so that the full expression is non-singular. Nevertheless, 
one might expect a simpler form with fewer such cancellations. Some such cancellations 
have already been carried out to arrive at cq. (6.7); besides the cancellation involving the 
unusual factor (6.6), individual recursive and overlap diagrams also contained factors of 
1/(5~| (3 + 4) |2~) , as in eq. (6.4). These squared factors were eliminated using spinor- 
product identities to combine and rearrange terms. 

Removing all of the factors of the form {a~\ {h-\- c)\d~) is unhkely to be desirable; after 
all, the simple forms of tree amplitudes found in refs. [9, 19, 20, 41, 42] are simpler than 
previous forms precisely because of the presence of some such denominators. However, the 
denominators present in eq. (6.7) occur in a rather asymmetric fashion, preventing the flip 
symmetry of the expression — symmetry under 

1^2, 3^6, 4^5 (6.8) 
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FIG. 7: The recursive diagrams for the amphtude ^7.1(1 ,2 , S"*", 4"*", S"*", G"*", T"*"), using the (1,2) 
shift. Diagrams (a), (b), (d) and (e) vanish. 




FIG. 8: The overlap diagrams for the seven-point amphtude. 

— from being manifest. Indeed, the factor (5~| (3 + 4) \2~) occurs in the denominator many 
times, yet its image under the flip (6.8), (4~|(5 + 6)|1~) never appears. 

We have carried out several numerical checks of cq. (4.28), after inserting into it the value 
of Rq from eq. (6.7). Besides verifying the absence of any spurious singularities, we checked 
that the amplitude satisfies the non-manifest fiip symmetry (6.8). We checked the multi- 
particle factorization in the S123 and S234 channels. (The S612 channel is related trivially to 
the S123 channel by the flip symmetry.) We also conflrmed the proper coUinear behavior, for 
real momenta, in all the independent channels S12, S23, S34 and S45. These checks leave little 
doubt that cq. (6.7) is the correct expression. 

Using our expression for the six-point amplitude as well as the seven-point cut terms from 
ref. [4] as input, it is then straightforward to evaluate the complete seven-point scalar loop 
amplitude, aS(1~, 2", 3+, 4+, 5+, 6+, 7+). The needed recursive diagrams are displayed in 
flg. 7, while the overlap diagrams are given in flg. 8. We have evaluated these diagrams and 
have conflrmed that the resulting expression for the amplitude has all the proper factorization 
properties in real momenta, and that all spurious singularities cancel. Combining this result 
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with the known A/" = 4 [3] and A/" = 1 [4] supersymmetric amphtudes, yields a complete 
solution for the seven-gluon QCD amplitude with the same helicity configuration. Although 
its construction is entirely straightforward, and parallels the six-point case, the seven-point 
result is rather lengthy, so we refrain from presenting it here. Since the original version of 
this paper appeared, the result has been presented [62], as a member of the infinite series 
of n-point amplitudes A^!1^{1~, 2~, 3+, 4+, . . . , n^), constructed using the methods of the 
present paper. 

These examples demonstrate the power of the factorization-bootstrap approach, system- 
atized here, as a complement to the unitarity-based method, for evaluating complete QCD 
amplitudes, including purely rational parts. The required diagrams are surprisingly simple 
to evaluate, not really more involved than tree-level diagrams. It is striking that what had 
previously been the most difficult part of a one-loop QCD calculation has been reduced to 
a simple computation. 

VII. CONCLUSIONS 

In this paper we presented a new method for computing the rational functions in non- 
supersymmetric gauge theory loop amplitudes. The unitarity method [3, 4, 5, 6, 7] has 
already proven itself to be an effective means for obtaining the cut-containing terms in 
amplitutes, so we may rely on this approach for obtaining such terms. To obtain the ra- 
tional terms we took a recursive approach, systematizing an earlier unitarity-factorization 
bootstrap [2]. 

Our systematic loop- level recursion uses the proof of tree-level on-shell recursion relations 
by Britto, Cachazo, Feng and Witten [10] as a starting point. There are, however, a number 
of issues and subtleties that arise, which are not present at tree level. The most obvious 
issue is that the tree-level proof relies on the amplitudes having only simple poles and no 
branch cuts; loop amplitudes in general contain branch cuts. Furthermore, as we have 
already discussed in refs. [11, 12], there are subtleties resulting from the differences of one- 
loop factorizations in complex momenta as compared to those in real momenta. These 
differences have important effects, unlike the tree-level case. At loop-level there are also 
spurious poles present, which would interfere with a naive recursion on the rational terms. 
In this paper we showed how to overcome these potential difficulties. 
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As an illustrative example of our approach, we described in some detail a computation 
of the rational terms appearing in the five-gluon QCD amplitudes with nearest-neighboring 
negative helicities in the color ordering, reproducing the results [47] of the string-based 
calculation of the same amplitudes. Although we did not describe it in any detail, we also 
confirmed that our new approach properly reproduces the other independent color-ordered 
five-gluon helicity amplitude. 

Next we computed the six- and seven-point QCD amplitudes 2", 3+, 4+, 5+, 6+) 

and 2~, 3+, 4+, 5"*", 6"*", 7+). The rational terms of these amplitudes had not been 

computed previously. Our computations of these terms use as input lower-point ampli- 
tudes [11, 47, 57, 58, 59], and the cut-containing terms of the amplitudes under considera- 
tion, obtained previously via the unitarity method [4] . For the six-point case we presented 
a compact expression for the complete amplitude. 

Another possible approach to obtaining complete loop amplitudes is via £)-dimensional 
unitarity [5, 13, 14, 17]. It would be worthwhile to corroborate the results of this paper 
starting from the known dimensional tree amplitudes [5, 13, 27]. It would be also be 
desirable to develop a first-principles understanding of loop-level factorization with complex 
momenta, instead of the heuristic one of refs. [11, 12]. 

The computation of rational function terms has been a bottleneck for calculating one- loop 
amplitudes in non-supcrsymmetric gauge theories with six or more external particles. We 
expect the technique discussed in this paper to apply to all one-loop multi-parton amplitudes 
in QCD with massless quarks. It should also work, without modification, for amplitudes 
that contain external massive vector bosons, or Higgs bosons (in the limit of a large top- 
quark mass), in addition to massless partons. Finally, we expect suitable modifications of 
the method to be applicable to processes with massive particles propagating in the loop. 
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